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Arguments are given that support the interpretation of the lattice QCD glueball and gluelump 
spectra in terms of bound states of massless constituent gluons with helicity-1. In this scheme, 
the mass hierarchy of the currently known gluelumps and glueballs is mainly due to the number of 
constituent gluons and can be understood within a simple flux tube model. It is also argued that 
the lattice QCD +_ glueball should be seen as a four-gluon bound state. 
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O ' I- INTRODUCTION 

o : 

QCD allows the existence of purely gluonic bound states, called glueballs, whose structure and properties deserve 
much interest (see Ref. [l[ for a recent review). An important theoretical achievement in this field has been the 
computation of the glueball spectrum in lattice QCD @, y, 0|- A plot of this spectrum is given in Fig. [TJ where the 
masses are expressed in units of the lattice energy scale r^ 1 in order to avoid larger error bars due to the determination 
qq ■ of this energy scale. Notice that r^ 2 can be understood as the energy density of a fundamental flux tube, that is, the 
' flux tube linking a quark-antiquark pair. Apart from lattice QCD, the glueball spectrum has also been computed by 
using effective approaches like Coulomb gauge QCD @ and potential models (see for example Refs. @, 0, H, M, HH ) • In 
potential models, glueballs are assumed to be bound states of constituent gluons. Within this framework, gluons are 
supposed to be massless or not, with either a helicity-1 or a spin-1. The physical properties of constituent gluons are 
actually still a matter of controversy, and even the constituent gluon picture could be questioned since potential models 
Qh| have serious difficulties in reproducing all the currently known lattice QCD data. It has been shown in Ref. [ll[ that 
the lightest glueballs, which have C = +, can be successfully modeled by a two-gluon system in which the constituent 
gluons are massless helicity-1 particles. The proper inclusion of the helicity degrees of freedom dramatically improves 
the compatibility between lattice QCD and potential models [ll|. This picture is also supported by the Coulomb 
<N ; gauge study of Ref. @. 

Motivated by these previous results, we propose in the present paper an interpretation of the whole glueball spectrum 
from a helicity-1 constituent gluon point of view. The relevance of this picture and its links to QCD are discussed in 
Sec. HH and the general principles of our classification scheme arc exposed in Sec. IIIII Glueballs and gluelumps are 
then analyzed in Sees . IIVI and fVl respectively. Gluelumps are actually bound states of the gluonic field in a static color 
octet source that have been studied first in lattice QCD [12 ]. Only a relative mass spectrum was computed in this last 
reference, while a determination of the absolute mass spectrum has been given in Ref. [l3| . The lowest-lying gluelump 
qq ' states are also plotted in Fig. [T] What we argue is that the low-lying gluelumps can be interpreted as bound states 
of a single constituent gluon. A simple flux tube potential model is then proposed in Sec. IVII in order to understand 
the mass hierarchy of gluelumps and glueballs. By combining the flux tube model and the helicity formalism, we will 
. [ finally compute in Sec. IVIII the mass of the lowest-lying H four-gluon state and compare it to the lattice QCD data. 
Conclusions will then be drawn in Sec. IVIlfl Appendix lAl gives some two-body helicity states not built in Ref. [ill ]. 
H ' while Appendix iBl recalls the way of computing the C-parity of few-gluon glueballs. 
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FIG. 1: Plot of the lattice QCD glueball spectrum (triangles). Data are taken from Ref. [If. The four lowest-lying gluelumps 
are also plotted (squares). Gluelump masses are taken from Ref. [l3| ]. Masses are given in units of the lattice energy scale r^ 1 . 

II. CONSTITUENT GLUONS AND QCD 
A. Effective approaches 

The basic assumption underlying our work is that a description of gluelumps and glueballs in terms of states with 
a given number of constituent gluons is relevant, or at least that it can be a satisfactory approximation of what these 
gluonic hadrons exactly are. The idea that low-energy QCD allows for an effective description of hadrons as bound 
states of constituent particles is not new: The classification of baryons and mesons with the quark hypothesis is a 
first historical example of the viability of such a picture. However, there should exist a way to connect the constituent 
approaches to QCD. To our knowledge, this longstanding problem is far from being completely solved, but interesting 
results have been obtained within different frameworks. We recall some of these approaches and comments on them 
in the following. 

First, we turn our attention to large- N c SU(7V c )-gauge theories, where N c is the number of colors. The 1/N C 
expansion in QCD, valid in particular at low energy, has revealed itself to be a powerful tool which has played an 
important role in numerous theoretical aspects of QCD (see e.g. Ref. for useful references). In the limit where the 
number of colors becomes very large, it has been shown that the description of baryons, made of N c valence quarks, 
becomes simpler (l5| . Model independent mass formulae can be obtained within this framework and lead to a very 
accurate description of experimentally known baryons, either light or heavy [l4| ; thus the large- N c limit appears to 
capture the essential features of QCD, which is characterized by "only" N c — 3. Moreover, it has recently been shown 
that the mass formulae obtained from a constituent quar k model are in remarkable agreement with those arising 
from the large- iV c expansion in the light baryon sector [161 ]. This comforts the constituent quark picture in the most 
difficult cases, when light (even massless) quarks are present. Of importance for our purpose is the recent observation 
that the number of valence gluons is a good quantum number for glueballs in the large- N c limit [13] • By analogy 
with the baryons, this suggests that the Fock space expansion of glueballs at N c — 3 is dominated by a particular 
component, characterized by its number of constituent gluons. In this spirit, the present work could be seen as an 
attempt to identify these dominant components for each known glueball and gluelump. 

Second, we can get useful informations from Coulomb gauge QCD [18j. In this approach, the QCD Hamiltonian 
is written with the gluonic field in the Coulomb gauge. A key feature of this gauge is that the elimination of the 
nondynamical degrees of freedom creates an instantaneous nonperturbative interaction. Properties of low-energy 
hadronic states are eventually computed from this Hamiltonian and the Fock space in which these states are defined 
is expanded on a quasiparticle basis, either quarks or gluons. A general property of hadronic states is then that their 
Fock space expansion converges quickly. This leads to a connection with constituent models since, again, it seems 
relevant to assume that a fixed given number of valence quarks or gluons can be associated with a given hadron. The 
notion of an instantaneous confining potential also naturally emerges from Coulomb gauge QCD. The connections 
between constituent quark models and Coulomb gauge QCD have been discussed in Ref. [19( while glueballs are studied 



in Refs. [20(. It is worth mentioning that in these last works, the low- lying glueballs are successfully described in terms 
of two-transverse gluon bound states, with a negligible three-gluon component. These conclusions are shared by the 
recent results of Ref. [III ], where pure two-transverse gluon bound states have been shown to accurately reproduce 
the lattice QCD and Coulomb gauge QCD data. 

In Coulomb gauge QCD as in many potential models, the constituent particles have a particular status since they 
are confined. In particular they gain a constituent mass, different from their bare mass, because of the confining 
interaction. One can look in particular at the Coulomb gauge study of Ref. |2ll | , where it is shown that constituent 
gluons, although having a zero bare mass, gain a running constituent mass which is about 0.7 GeV at zero momentum. 
They become then distinguishable from the exchanged gluons, which mediate the interaction without being themselves 
confined. 

An effective QCD Hamiltonian for various hadronic systems can also be derived from the field correlator method 
f22| . In these models, the picture of constituent particles (quarks or gluons) linked by color flux tubes emerges 
naturally. For instance, the constituent mass of the confined gluons is defined as (p 2 ), that is their average kinetic 
energy. Interestingly, this quantity is also typically of about 0.7 GeV for the glueball ground state. 

B. The correlator method 

We have just shown that several effective approaches of low-energy QCD agree with the constituent gluon picture 
for glueballs. Even if the Fock space expansion of a glueball or gluelump may be nontrivial, these models actually 
suggest that one component with a given number of valence gluons may be dominant. Now we make some comments 
on the correlator method, that allows to compute glueball masses in lattice QCD for example |23j . 

Let us first consider the initial state (at zero euclidian time) 1^(0)) = $^(0) |0), where ^^(O) is an operator 
creating a glueball or a gluelump out of the QCD vacuum |0), and where R denotes the irrep corresponding to particular 
J quantum numbers. Although highly nontrivial, the evolution of this state is computable, with lattice methods 
in particular [H, H3|- After a time r one gets the state |*(r)) = <E> (fl) (r) |0), and the correlator C(r) = ( v E'(t)| <3>(0)) 
can be evaluated. As intuitively expected, lim^oo C(r) cx e~ ma T , rog being the mass of the lightest state which can 
be created by the operator $(- R ) . Examples of such operators are given in Table [IJ A rule giving the mass hierarchy of 
glueballs has been formulated in Ref. [2a|: The higher the mass dimension of an operator is, the heavier the resulting 
hadron is. For example, the lowest I H gluelump operator has dimension 2, while the first ++ and 1 H glueballs have 
dimension 4 and 6 respectively. Figure Q] shows that the observed lattice QCD mass hierarchy follows this dimension 
rule. There are two ways of increasing the dimension of an operator: Either to increase the order of the operator in 
the chromo-electric or -magnetic fields E a and B a respectively, or to add a covariant derivative D l (as in the 2~~ 
gluelump) . Both ways have a very natural explanation within a constituent gluon picture. Indeed, a transverse gluon 
creator is contained in each E a or B a factor. Assuming thus that each such factor generates a constituent gluon, we 
find that gluclumps are made of a single constituent gluon, and that they are logically lighter than glueballs, which 
are made of at least two constituent gluons. Another illustration of this interpretation is that a C-odd glueball can 
be created with at least a cubic expression in the fields: As possible three-gluon states, they logically appear to be 
heavier than two-gluon states. Within this picture, the addition of totally symmetrized covariant derivatives does not 
affect the number of constituent gluons, but adds momentum: The covariant derivative rather generates excitations 
of the considered few-gluon system (25j . 

Even if ^(O)) can be seen as a state with a given number of constituent gluons, because it is created by operators 
such as those of Table HI the final physical state \^(t — » 00)} could have a far more complicated Fock space expansion 
due to the complex QCD interactions that come into play. But, as mentioned above, this does not seem to be the case. 
At this stage a connection with potential models can be made: They all rely on the assumption that it is possible 
to find an appropriate effective potential satisfactory reproducing the properties of physical hadrons (« |\P(r — » 00))) 
seen as systems of constituent par ticles (« |^(0))) interacting via this effective potential. To conclude, it is worth 
mentioning the results of Ref. [26|, where it has been shown that the lattice QCD ++ glueball wave function is 
compatible with the wave function of a flux tube Hamiltonian of the form 1\ffi 2 ~ + ar — n/r in the ++ channel. 
We recover the idea that the scalar glueball can be seen as a dominantly two-gluon state, where the two constituent 
gluons (« 2^JjP) interact via the effective potential ar — n/r, 

III. GENERAL PRINCIPLES 

Following the above discussion, it seems very reasonable to build our study from this assumption: Glueballs and 
gluelumps can be seen, at least in first approximation, as bound states of constituent gluons, these being massless 



TABLE I: Examples of operators creating glueballs and gluelumps with a given J . The () denote a complete symmetrization. 
The repeated color indices a, b, c are intended to be summed and the superscripts i, j are space indices. 
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particles with helicity-1. The quantum states describing bound states of few-gluon systems can be efficiently formulated 
within the helicity formalism, introduced in Rcf. [27j . Hereafter, we recall the main points of this formalism. Let us 
denote \tp(p, A)) the quantum state of a particle with momentum p and helicity A. If this particle is massive of spin s, 
then A can take the 2s + 1 allowed values for a spin projection along a given axis. But, if the particle is massless, only 
A = ±s is allowed. The question is now: How to write a two-particle helicity state in the rest frame of the system 
such that both particles are coupled to a given value of the total spin J with a projection J z — M and with a given 
parity? The general answer is given in Ref. [27| . from which it can be deduced that the helicity states 

|A!,A 2 ; J p ,M,e) = -± {Sl J MM _ X2 [|V(?W> ® l<K-£A 2 )>] 

+^j IM _ Xl \\m -AO) ® w-a - A 2)}]} (i) 
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d0 sin P^, A (0, 0, -0) 0, -0) 0), (2) 



have the desired features, with R{a, [3, 7) the rotation operator of Euler angles {a, (3, 7} and V J M x (a, (3, 7) the Wigner 
D-matrices. Let us remark that | Ai, A 2 ; J p , M, e) = e | — Ai, — A 2 ; J^, M, e); both states are equivalent up to an 
irrelevant e factor. 

The parity of the state ([T]) is given by P = eriiri2(—l) J ~ Sl ~ S2 T]i and Si being the intrinsic parity and spin 
of particle i respectively. The parameter e, which can be equal to ±1, fixes the final value of the parity. It is worth 
mentioning that, following the usual rules of spin coupling, 

J>|Ax-A 2 |. (3) 

Another important point to notice is that when both particles have a spin degree of freedom, the helicity basis, 
spanned by the helicity states ([1]), is equivalent to an usual | 2S+1 Lj) basis j28|. When at least one of the particles 
is massless, both basis are no longer equivalent but the helicity states can still be expressed as particular linear 
combinations of | 2S+1 Lj) states [27j . In order for our discussion to be complete, we have to mention that, in the case 
where the two constituent particles are identical, the helicity states have to be totally (anti)symmetrized, following 
the bosonic or fermionic nature of the particles. The action of the permutation operator on the helicity states can be 
found in Refs. [HI, and we will not mention it here explicitly. But, it is worth saying that the symmetrization of 
the systems leads to constraints on the total angular momentum, thus to particular selection rules. This fact will be 
illustrated in the case of two-gluon glueballs. 

A nice property of the quantum state {T]) is that it is well defined for massless constituent particles. In particular, 
for a system made of two gluons, we have = ±1 and P = e(— 1) J . An important quantity that can be computed in 



this case is the average square orbital angular momentum, that reads (see Ref. ll[ and Appendix [A"]) 



L 2 ) = J(J +l) + 2AxA 2 . (4) 



A first idea to classify the glueballs is that \L 2 J roughly sets the energy scale of a given state: The more rotational 

energy is contained in a glueball, the more the state is heavy. Notice that Eq. (0]) is valid for two-gluon states only. 

The prescription to write a three-body helicity state has firstly been given in Ref. [29| and is based on an intuitive 
recoupling scheme. The first step is to write the two-body helicity state corresponding to, say, particles 2 and 3 in 
the center of mass frame of these two particles. This is done thanks to the result ([1]). Then, the Lorentz-boosted 
helicity state of the (2,3) cluster can be coupled to particle 1 with the same technique. The boost is needed because 



one wants the three-body helicity state to be expressed in the three-body rest frame, not in the rest frame of the (2,3) 
cluster. We propose the following notation to denote a three-body helicity state 

|(A 2 ,A 3 ;r,"i,e),Ai;J p ,Af,p). (5) 

This last ket has to be understood as follows. Two constituent particles with helicities A2 and A3 are coupled to a 
given value of j > |A 2 — A3 1 , with a parity it — e?72?73(~l) : ' S2 ~ S3 and a helicity m following Eq. ([I]). Then the (2,3) 
cluster, seen as a pointlike massive particle with spin j and helicity m, and the last gluon with helicity Ai are coupled 
to the total value of J > \m — Ai|, with a parity P = prii-K{~l) J ~i~ Sl — ep?7i?72'73(~l)' 7 " 1 " 2 " 3 - This time, the 
product ep determines the total parity, with p = ±1. For a three-gluon system, the parity is given by P = ep(— l)' 7 . 
Although formally clear, the situation is much more complicated for three-body systems, because an infinity of helicity 
states of the form ([5]) can theoretically lead to a given J p state. Again, a classification can be made by assuming that 
a J PC three-gluon glueball is dominated by the helicity state for which the value of both the square orbital angular 
momentum of the cluster and the total square orbital angular momentum are minimal. 

We propose to adopt a similar procedure to study four-gluon states: First, to write the helicity states for the (1,2) 
and (3,4) clusters and, second, to couple these two clusters. The advantage of doing so is that both clusters must 
be coupled as massive particles. Consequently, the usual | 2S+1 Lj) basis can be used in this last step, and the effect 
of helicity is restricted to the two-gluon helicity states, which have been extensively studied in Ref. [ll|. Such a 
procedure can actually be extended to many-gluon or many-body systems, with an obvious increase of complexity 
with the number of constituent particles. Moreover, a global symmetry can become very difficult to implement when 
more than two constituent particles are identical (see for example Ref. (30j , where the case of three identical fermions 
is considered). 

Group theory generally gives powerful tools concerning the classification of hadrons and it will also be interesting 
to use it in the following. One could represent a constituent gluon by the Young diagram a , the box denoting the 
fundamental representation of the little group SO (3) [3l| and a = 1, . . . , n being the color index associated with the 
gluon. Obviously n — — 1 with N c the number of colors. The isospin of a gluon is zero; this degree of freedom is 
thus trivial and will not be explicitly taken into account here. We then expect that the lightest glueballs made of N g 
constituent gluons will be those for which the J PC quantum numbers can be generated by the intrinsic color and spin 
degrees of freedom, i.e., since gluons are bosons, by the states appearing in the decompositions of the symmetrized 
tensor product 

(□ Hl ® ■■■® n^"): (6) 

where the symbols () denote the total symmetrization. Only the total spin J is a relevant observable, but the helicity 
states can always be decomposed in a usual | 2S+1 Lj) basis. The J PC states appearing in the decomposition Q will 
consequently be those for which a S'-wave component is present in the helicity state. Thus, we can expect that this 
will correspond to the lightest glueball states with N g gluons. Before illustrating these general ideas by considering 
explicit examples, we refer the reader to Ref. [321 ] for a discussion about representation theory. More precisely, the 
symmetrized tensor products - demanded by the Pauli principle - that arise when doing tensor products of several 
identical boxes Q a are referred to as plethysms. It is worth mentioning Ref. [Hj], in which a fully group-theoretical 
classification of glueballs is performed following a SU(3) x SU(2) scheme. But, the influence of gluon helicity through 
the helicity formalism has, to our knowledge, never been taken into account and leads to interesting and new results 
as we will see through this paper. 

IV. THE GLUEBALL SPECTRUM 
A. Two-gluon states 

Bound states of two gluons are the most studied purely gluonic systems in the literature. Since their properties 
are rather well-known, it is useful to re-analyze them using the method introduced in Sec. IIIII in order to check its 
compatibility with previous works. Let us begin by a Young diagrams analysis. For a system made of two constituent 
gluons, one has 

(□ a ®n fc ) = {m©«} (ab) ©n [ahl , (7) 

where (a±a2 ■ ■ ■ o-N g ) and \a\a2 ■ ■ ■ ttjvJ denote respectively the totally symmetric and antisymmetric color wave func- 
tions. The notation {} allows avoiding to repeat the same color superscript for different spin irreps, and the • stands 
for a spin zero state. Let us briefly make a comment about the dimension of the symbols appearing in this last 



equality. The Young tableau Q a has a dimension 3n, i.e. the product of the spin and color degrees of freedom. Then, 
the left hand side of Eq. ([7]) has the dimension 3n(3n + l)/2 since only its totally symmetric part is considered. The 
three terms appearing on the right hand side have dimension 5n(n + l)/2, n(n + l)/2, and 3n(n — l)/2 respectively, 
the sum of which gives 3n(3n + l)/2 as expected. 

Since a two-gluon state must be in a color singlet, its color wave function is trivially [8,8] la , the resulting singlet 
being symmetric under the permutation of both gluons. Only the first two terms of the decomposition ([7J then 
correspond to colorless glueballs: An antisymmetric color wave function cannot correspond to a color singlet. Notice 
that the rj[ Qb l state could correspond to colored glueballs with the color wave function [8,8] 8a , [8,8] 10a or [8, 8] 10a , 
and could become relevant in a deconfined medium as the quark-gluon plasma. Notice that the superscript X a / S is 
used to denote a symmetric (s) or antisymmetric (a) representation X. The first two states of the decomposition ([7]) 
are J c = 2 + and + glueballs respectively. The C-parity of a two-gluon system is always positive as shown in 
Appendix [5] We do not discuss the parity of these state at this stage since only the helicity formalism can bring 
relevant information about it. 

The helicity states corresponding to two-massless gluon systems can be explicitly written thanks to Eq. (fT|). Since 
all the necessary calculations have been performed in Ref . , we only mention here the final results concerning the 
two-gluon helicity states and their properties. We have 

|l,l;(2fc)+,M,l) =► 0++, 2++, 4++,... 

|l,l;(2fc)-,M,-l) => CT+, 2"+ 4"+, . . . (8) 
|l,-l;(2fc + 2)+,M, 1) 2++, 4++, 6++,... 

|l,-l;(2fc + 3)+,M,-l) => 3++, 5++, 7++,..., 

with k — 0, 1, ... The selection rules that appear on the total spin J are actually a consequence of the fact that these 
helicity states have to be totally symmetric [ll|. In particular, the state |1, — 1; 2~, M, — 1) is forbidden for such a 
reason. In agreement with Eq. (J7J and Fig. [TJ the lightest glueballs correspond to |1, 1; + , 0, 1} and |1, —1; 2 + , M, 1} 
respectively, for which the square orbital angular momentum is 2 and 4 respectively. The ^ glueball, given by 
1, 1; CP, 0, —1), is heavier than the ++ one although its square orbital momentum is also equal to 2. This shows that 
other effects come into play into glueballs: We suggested in Ref. [ll[ that the mass splitting between the scalar and 
pseudoscalar glueball is due to instanton- induced interactions. Finally, always by looking at Eq. we are led to the 
conclusion that the 2 h is lighter than the 3 ++ one, as observed in lattice QCD (see Fig. [T]). 

An important remark is that no J PC = 1 ±+ glueball should be present in the energy range of the currently observed 
C = + states because of the selection rules on the states ©. This is indeed the case in lattice QCD and gives support 
to the framework that we propose. We notice that this result is obviously similar to Yang's theorem stating that a 
vector meson cannot decay into two photons [341 ] - But, if constituent gluons are spin-1 particles, low-lying 1 ±+ states 
are then allowed, in disagreement with what is found in lattice QCD [ljj - consider the \ 3 Pi) state for example. The 
presence of such vector glueballs is actually due to the fact that Xi = is allowed with spin-1 constituent gluons. 



B. Three-gluon states 



When three constituent gluons are taken into account, the following decomposition holds 

(□ a ® n b ® lt) = {nnn © u} (abc) © {m © ni^ 1 ' © • [abc \ (9) 

12 11 

where M b ' denotes a color configuration with the mixed symmetry 

can be made with this mixed symmetry (see Appendix IB"]). So, the third and fourth states of the decomposition @ 
cannot be colorless three-gluon states. But, the symmetrical color wave function (abc) corresponds to the following 
color singlet, [[8,8] s ,8] ls , and leads to C = — glueballs. Similarly, the antisymmetrical color wave function [abc] 
corresponds to the following color singlet, [[8,8] 8a ,8] la , and leads to C = + glueballs as already pointed out in 
Ref. [2a]. Consequently, the glueballs corresponding to the first two states of decomposition are J c = 3~ and 1~ 
states respectively, while the last state is a + one. We will actually focus on the C = — three-gluon states in this 
section since no mixing of these states with two-gluon ones can exist. We remark that C = — three-gluon glueballs 
are formally identical to three-photon states. 

Following the arguments that we gave in the previous section, the lightest three-gluon states should be those for 
which the two square orbital momenta appearing in the three-body system have both a minimal value. Ignoring 
possible instanton effects, the most obvious possibility to build such a system is the following 

5|(l,l;0 e ,0, e ),l;l-^,M,p), (10) 




. But, it can be checked that no color singlet 



where j and J have the minimal allowed value, leading to 1 ±_ glueballs, and where S is the three-body symmetrizer. 
We stress that, formally, an infinity of {j, m, e, p} sets can lead to a given J PC state. For example, the quantum state 
of the l" 1 glueball is rigorously given by a linear combination of all the totally symmetric three-body helicity states 
for which the values of {j, m, e, p) can lead to J PC = l" 1 state. Here we assume that the lowest-lying glueballs are 
dominated by the component in which the minimal amount of rotational energy is present. Similarly, light J PC = 3 ±_ 
glueballs can a priori be built from the following state, 

S\(l,-l;2+,-2,l),l;3-P,M,p) , (11) 

where again both j and J are minimal following the values of and m that have been chosen, and where the rotational 
energy is the lowest possible for the cluster (2,3). 

In order to understand the mass hierarchy of the lattice QCD states, we propose to introduce the notion of natural 
and unnatural parity for glueballs. It is indeed known in hadronic and nuclear physics that states with natural parity 
are lighter than those with unnatural parity. We suggest that the natural parity of a glueball made of N g gluons is 
given by (— l)^ 9 (— 1) J , where the first factor comes from the intrinsic parity of the N g gluons, and where the second 
one is the expected parity of a state with spin J. For two-gluon glueballs, the ++ and 2 ++ states have a natural 
parity, and are lighter than the h and 2 h ones. In the case of three-gluon systems, the introduction of natural 
parity leads to the prediction that the 1 H (3 H ) glueball should be lighter than the 1 (3 ) one. Thus, we predict 
the following inequalities concerning the glueball masses 

M 1+ -<Mx—, M 3 +-<M 3 — . (12) 

This is indeed what is found in lattice QCD, as one can see in Fig. [T] We mention also the inequality M\±- < M 3 ±-, 
which can be deduced from the angular momentum content of the states (|TU)l and (jTTJ) . 

An important point to notice is that no low-lying ±_ glueball can be obtained since the lowest allowed value for J 
is 1. This result has already been mentioned in Ref. (36[, where it is shown that no (pseudo)scalar three-photon state 
can exist, excepted in very asymmetric configurations that seem irrelevant to describe bound states. The 0" 1 glueball 
that is observed in lattice QCD should then be seen as either a highly excited three-gluon state or as a low-lying 
state made of at least four-gluons (two gluons are forbidden by the negative C-parity). In both cases its large mass 
is explained (it is the heaviest glueball that is currently known in lattice QCD). But, as we will see in the following 
sections, the hypothesis of a low-lying four-gluon state appears to be particularly relevant. We also remark that, 
when spin-1 constituent gluons are assumed, the trivial state |(0, 0; £ , 0, e), 0; 0, e) , can lead to a glueball 

which should be among the lightest ones in the PC = channel, but which is not seen in lattice QCD. This favors 

the helicity-1 approach since the C = — sector is known in a quite large mass range, and no glueball has been 
observed yet. It is finally worth mentioning that with spin-1 gluons, a H three-gluon glueball is also allowed, but it 
has been shown in a previous work that its mass is far too close of the other three-gluon candidates with respect to 
the mass splittings that are observed in lattice QCD [lj|. 

The additional 2 and 2 H states that are present in Fig. [1] could be understood as the first excitations of the 
helicity state (|10[) for example. By excited state, we mean a state whose total spin J is higher than the lowest allowed 
value. Since the parity is proportional to (— 1) J , the 2 (2 H ) glueball could be an excitation of the l" 1 (1 ) 
glueball. This point is still compatible with lattice QCD since the 2 (2 H ) glueball is observed to be heavier 
than the 1 H (1 ) glueball. Notice that M 2 -- < M 2 +- is observed: The glueball with natural parity is still the 
lightest. It could be argued that states such as |(1, —1; 2 + , 1, 1), — 1; 2 P , M, p) could lead to 2 ±_ glueballs without any 
excitation. But the problem actually comes from the value m — 1 of the cluster helicity. Indeed, if m — ±1, low-l yin g 
glueballs such as |(1, —1; 2 + , 1, 1), 1; P , 0, p) are allowed. This would be in contradiction with the results of Ref. [3a |. 
forbidding (pseudo)scalar totally symmetric three-gluon systems, or at least low-lying ones. We actually think that 
the case \m\ = 1 is forbidden by the requirement that the three-body helicity states have to be totally symmetrized. 
A confirmation of this point needs an explicit implementation of the three-bod y sy mmetrizer and a computation of 
its action on the helicity states. Such a task is far from being trivial (see Ref. [3fJ for a fermionic example), and is 
out of the scope of the present paper. Nevertheless, we suggest that 2 ±_ three-gluon glueballs are most likely excited 
states. 



C. Four-gluon states 

The situation becomes a little bit more involved in the case of four-gluon systems. One obtains 

(□ a ® n\ b ® lt ® n d ) = ( i i i i i © m © •} {abcd) © {m © 

© □ M -- 11 © { rm © m © □ } M ^ , (13) 



and 



respec- 



where M^^j, Af^bcd and Af^^j 1 ' denote the color configurations with the symmetries | | I I 

tively. It is worth discussing a bit the different allowed color singlets in the case of four-gluon systems. By using a 
(1,2)— (3,4) recoupling scheme as it is done within the helicity formalism, we can write all the confined color singlets 
as follows: [[8, 8] 8 «, [8, 8] 8 *] 1 , [[8, 8] 8 «, [8, 8] 8 -] 1 , [[8, 8] 8 % [8, 8] 8 -] 1 , [[8, 8] 15 - , [8, 8] 10 -] 1 , and [[8, 8] 27 = , [8, 8] 27 *] 1 (the 
[[8, 8] 1 , [8, 8] 1 ] 1 color wave function would lead to a deconfined four-gluon state which would obviously split into 
two two-gluon glueballs). These color singlets lead to different values of the C-parity that have been computed in 
Appendix [Bj We remark first that no colorless glueball can be made with the M' 3 ' 1 ' symmetry. Only the first six 
irreps of the decomposition (|13p are thus relevant to our purpose. The corresponding states can be identified with 
J c = {4, 2, 0} + , {2, O}*, and 1~ glueballs respectively. Due to the many configurations that are possible, an exhaus- 
tive analysis of four-gluon states would be rather tedious, and not very worthy since only the 0" 1 glueball appears 
as a relevant four-gluon candidate in lattice QCD. In this section, we rather show that a four-gluon H glueball is 
allowed and we focus on its specific properties. 

It is known that, at tree level, the short-range interactions between two constituent gluons are mainly proportional 
to C gg — 6, where C gg is the quadratic Casimir of su(3) in the representation of the gluon pair [37| • The configurations 
in which a gluon pair is in the 10 a or 10 a representation would lead to vanishing short-range interactions of one- 
gluon-exchange type, while the 27 s representation would lead to repulsive short-range forces. Color wave functions 
in which the gluon pairs are in a color octet seem thus favored since the short-range interactions are attractive in 
this case. Configurations involving the 27 s representation should conversely lead to heavier glueballs because of the 
existence of repulsive forces. It can be checked that, among the color wave functions involving color octets only, 
[[8, 8] 8a , [8, 8] 83 ] 1 is the unique configuration that has a negative C-parity. But, the resulting color singlet can only 
have the mixed permutation symmetry M' 2 ' 1 ' 1 ', that cannot be associated with a low-lying J = glueball as suggested 
by Eq. (fT3"]) . A light H glueball can consequently not be built from color configurations where each gluon pair is 
in a color octet; higher representations are needed. As it is shown in Appendix [B] a color wave function of the type 
[[8, 8] 10 °, [S^] 10 "] 1 can lead to a C = — glueball with the Afl 2 ' 2 ! permutation symmetry. Consequently, the unique 

| 2 2 | 

±_ state that can be built from Eq. (fl~3|) corresponds to • M <jf,'cd. Let us explicitly build a H state by using the 
helicity formalism. 

One has first to remember that the 0^ four-gluon glueball should have a particular mixed helicity symmetry of the 
Let us briefly come back to the color wave function. It reads [[8i, 82] 10a , [83, 84] 10a ], where indices labeling 



form 

the constituent gluons have been added. In this recoupling scheme, the helicity states must have precisely the mixed 
in order for the Pauli principle to be satisfied. Both gluon clusters are thus totally antisymmetric, 



symmetry 

a configuration that is forbidden for two-gluon glueballs. Let us define the antisymmetrizer of the cluster by 
A = (1 — Pij)/2, Pij being the permutation operator. Then it can be deduced from the results of Ref. [ll[ that the 
action of the antisymmetrizer on the two-gluon helicity states is given by 

2A\l,l;J p ,M,e) = [l - (-1) J ] |l, 1; J p , M,e) , (14a) 
2A\l,-l;J p ,M,e) = [l — e(— 1) J ] 1 1, — 1; J p , M, e) , (14b) 

leading to selection rules constraining the total angular momentum. Remember that we are looking for the lightest 
H state: Both gluons clusters should be as light as possible. By looking at Eq. Q, one can observe that the lightest 

antisymmetric two-gluon clusters correspond to the states |1, 1; l _e , M, e) and |1, — 1; 2~, M, —1), for which (^L 2 ^ = 4. 

Such a vector state is coherent with Eq. ([7}. 

Now that the quantum state of the clusters is known, we can lastly couple them in order to obtain a 0^ state 
(the negative C-parity is guaranteed by the color wave function). In order to satisfy at best the mixed symmetry of 
the four-gluon state, we will ask for the two clusters to be identical and for the two-cluster state to be the lightest 
J = one which is totally symmetric. This state is simply the [Sq) one, either spin-1 or spin-2 clusters are assumed. 
Notice that the parity of the two-cluster state is (±1) 2 (— 1) L = +1 as demanded. This is also the natural parity of a 
J = four-gluon glueball. 



V. A LINK WITH GLUELUMPS 



As it was recalled in the introduction, gluelumps are bound states of the gluonic field in a static color octet source. 
This color octet source is assumed to have the quantum numbers of the vacuum, i.e. J PC = ++ . Following the 
calculations of Ref. [l3T |. the lowest-lying gluelumps are lighter than all the currently known glueballs in lattice QCD. 
In consequence, we propose to identify the lowest-lying gluelumps with states made of only one constituent gluon 



(the static source holds for a scalar particle in a color octet, fixed at the center of mass). Such one-gluon states are 
generated by the tensor product 

□ ^.^□HffiQH, (15) 

where • represents the color octet static source. Note that no symmetrization must be achieved between a gluon 
and a static color source. As seen in Appendix [Bj only the symmetric color state (ab) can be reached by a color 
singlet, with C = — this time. The corresponding helicity states are given by 1 1, 0; J p , M, e), with a parity given by 
P = e(— 1) J as the static source is assumed to have a positive intrinsic parity. We notice that J > 1, so that the two 
lightest gluelumps should be the 1 H and 1 ones, followed by their first excitations, i.e. the 2 and 2 H states 
respectively. If the lattice QCD gluelump states are ordered by increasing mass, one finds the following quantum 
numbers: l" 1 , 1 , 2 , 2 H , 3 H , ++ , . . . 12]. This ordering is coherent with the one-gluon picture and with the 
natural parity criterion until a ++ state is reached, because such a positive C-parity state cannot be reached with a 
single constituent gluon. 

As pointed out in Ref. [38[, a ++ gluelump should at least contain two constituent gluons. Let us check this point 
within our formalism. The lowest-lying two-gluon gluelumps are expected to be generated by the following tensor 
product 

(□ a ©U fc )®- c = {[Jffi •} {ahc) © {m © □ © •} M ^ 1 ' © U labc] , (16) 

where the •( abc ) state with the natural parity should be the lowest-lying two-gluon gluelump. It indeed corresponds 
to a ++ state of the form 

5 23 |(l,l;0 e ,0, e );0;0+,0,e), (17) 

where the two gluons have been coupled first and symmetrized. Such a state has the lowest possible values for both 
the internal and total angular momenta and is similar to a two-gluon glueball. In particular, its mass is equal to 
roM ++ — 5.02 ± 0.46 [1J], that is a mass which lies in the typical range of the low- lying two-gluon glueballs (see 

Fig.m). 

When gluons have a spin-1, a light one-gluon gluelump can be simply built from |0, 0; 0~, 0, — 1), that is a 
single gluon with the zero spin projection. Interestingly, no gluelump has been observed yet in lattice QCD, still 
favoring the helicity- 1 picture that we develop in the present paper. 



VI. MASS HIERARCHY OF GLUEBALLS AND GLUELUMPS 



We have observed through the previous sections that the different lowest-lying one-, two-, three- and four-gluon 
states can be separated into glueballs with natural parity or unnatural parity. The corresponding low-lying glueballs 
that are currently observed in lattice QCD are listed in Table. [TTJ The left part is concerned with the natural parity 
glueballs (P = +) while the right part is concerned with the unnatural parity states (P = — ). In order to explain the 
mass hierarchy of these states, we define the spin-averaged mass as usual by 



TABLE II: Left: Natural parity low-lying iVg-gluon candidates. The quantum numbers and the masses of these states are 
denoted by J PC and tqMjpc respectively. Averaged masses are computed following Eq. (|18p . The N g = 1 states are taken 
from Ref. p]| while the other data are taken from Ref. Q. Masses are given in lattice units. Right: Idem for the unnatural 
parity candidates. 



Ng r° r MjPc [4, 13] r M(JV g ,l) t q Mjpc [4, 131 r M(N g , -1 



1 + - 
0++ 
2++ 
1+- 

3+- 

o+- 



2.25±0.39 
4.16±0.15 
5.83±0.11 
7.27±0.11 
8.79±0.12 
11.66±0.19 



2.25±0.39 
5.55±0.12 

8.33±0.12 

11.66±0.19 



1— 
0"+ 
2"+ 
1" 
3~ 



3.18±0.41 
6.25±0.12 
7.42±0.14 
9.34±0.13 
10.25±0.14 



3.18±0.41 
7.23±0.14 

9.98±0.14 
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FIG. 2: Plot of the spin-averaged masses of the natural parity (full circles) and unnatural parity (empty circles) states of 
Table |ll] versus the presumed number of constituent gluons. Equation (|19[) is plotted with the optimal values (|20[) for (3 and 7. 

where Mjpo is the mass of the glueball which is a Ag-gluon candidate with J PC quantum numbers. The case N g = 1 
is the gluelump case, while N g > 1 corresponds to glueballs. It appears that the spin-averaged masses grow linearly 
with the presumed number of constituent gluons; we then propose to fit the data of Table UT1 with the following form 

r M(N g ,P) = 0N g - 1 P. (19) 

A linear regression on the natural parity states leads to 

= 3.10 ± 0.08, and 7 = 0.81 ±0.21. (20) 

As it can be observed in Fig. [21 formula (|19p with the fitted values (|2"0")) reproduces rather well the spin-averaged 
masses. It is worth mentioning that this last mass formula predicts a four-gluon glueball in the mass range 
r M — = 13.21 ±0.53. 

Formula (TT^J) actually gives support to the constituent gluon picture, each constituent gluon roughly bringing to 
the total glueball mass an amount which is equal to 0. Additional mechanisms (instanton-induccd interactions, spin- 
dependent forces, etc.) are expected to come into play as higher-order terms and are all contained in the -7P term. 
If the particular value of 7 probably needs a very accurate model to be understood, can be obtained within the 
framework of the flux tube model. Indeed, can be simply obtained from formula (|19p 

= [roM(N g , 1) + r M(N g , -1)] . (21) 

The most simple case is the gluelump one, for which N g = 1. is then related to a typical gluelump mass, that can 
be obtained from a gluelump Hamiltonian in the flux tube model. At the dominant order, it reads 

H= \f]p+ar~?>— , (22) 
r 

where p and r are conjugate variables and where r — \f\ is the separation between the constituent gluon and the 
static source. It is spin-independent but, in our formalism, the effect of the gluons helicity is included at the level 
of the wave functions, not in the Hamiltonian itself. Although we do not use an explicitly covariant formalism, the 
semirelativistic kinetic term \pjp" should be relevant to describe the dynamics of the constituent gluon. a is the 
energy density of the flux tube linking this gluon to the static source, and a s is the strong coupling constant. The 
flux tube actually comes from the nonperturbative part of the QCD interactions and generates the confinement, while 
the Coulomb part is the dominant part of the one-gluon-exchange potential between the constituent gluon and the 
static source [33]. Assuming the Casimir scaling hypothesis, one can write a — (9/4)er with a the energy density of 



the flux tube in a meson. In order to compare our results with lattice QCD, we can define a — 1/r" 2 , and work with 
the new coordinates q = rgp and x = r/r$. Then, Hamiltonian (|22p becomes 

r H = . (23) 

4 x 

Following Eq. (f2"Tj) and Table HU (3 should be given by (M]+- + M 1 --)/2, both masses being computed from this 
last Hamiltonian. As mentioned in Sec. El two helicity states correspond to a single gluon with helicity-1: \^±) = 
|l,0; J p f M, ±), with J > 1 (see Appendix IA"]) . The eigenstates of Hamiltonian can be accurately computed 

using the Lagrange-mesh method [39[, using (^ili 2 |\I>±) = J(J + 1). The ground state of Hamiltonian (f2"3")) is 
then obtained for J = 1 and n = (n is the radial quantum number), and it is clear that M\+- = M x -- with 
this simple model. Let us denote M n j the eigenvalues of Hamiltonian (|23|) . One gets, after a numerical evaluation, 
Mot — 3.09 p» /3 for a s = 0.32. This value of the strong coupling constant is quite usual in potential models. This 
achieves to clarify the physical content of the coefficient (3. An important point to stress is that, in a usual, spin-based 

model, the lowest allowed value for (l 2 ^ is and would lead to the conclusion that (3 has to be identified with Afoo- 

But, it is readily checked that Mqq — > — oo unless a s = 0, which is not physically relevant. In this case, the use of 
spin-1 gluons clearly leads to unphysical behaviors. 

An estimation of (3 can also be obtained by considering the case N g = 2. Particularized to two-gluon glueballs, the 
flux tube Hamiltonian (|23[) simply becomes 

r # = 2^+ ?a - 3— , (24) 
4 x 

were the kinetic parts of two gluons have been included. In the same way than for gluelumps, f3 can be computed 
from the eigenstates of Hamiltonian l|24j) . In particular, the knowledge of the ±+ and 2 ±+ glueball masses finally 
lead us to a well-defined value for j3. We refer the reader to Ref. [ll[ for a detailed discussion about the computation 
of the two-gluon glueball spectrum within the helicity formalism. Once all calculations are done, with a s — 0.32 as 
for gluelumps, we find /? = 3.33, which is not far from the optimal value. 



VII. THE 0+ FOUR-GLUON GLUEBALL 



We have outlined in Sec. IIV CI the fact that the + glueball, which is currently the heaviest one that has been 
observed in lattice QCD, is a relevant four-gluon glueball candidate. We have also shown that such a system has 
a particular structure: Its color wave function should be [[8, 8] 10a , [8, 8] 1011 ] 1 - (see Appendix IB")) , with a mixed 
permutation symmetry of the form 



a lattice QCD glueball is proposed 

quite crude but physically relevant four-body calculation of the + 
the helicity formalism. 



To our knowledge, it is the first time that a four-gluon interpretation of 
A confirmation of this point is then important, in particular by performing a 

mass within a constituent gluon approach and 



A. Hamiltonian 



As in Hamiltonians used above, we adopt the flux tube picture stating that straight flux tubes start from the 
constituent particles and meet in such a way that the energy of the system is minimal. It was suggested in the 
recent lattice QCD studies of Refs. [40( that the less energetic configuration for many-gluon systems (more than 
two) is the one in which each constituent gluon is linked to its two nearest neighbors by two fundamental flux 
tubes. Phenomenological arguments also support this picture |4lj |. In agreement with this idea, we suggest that the 
Hamiltonian describing a glueball made of four gluons is given at the dominant order by 

4 

H = Y,vfi +a E i*-*<+ii, ( 25 ) 

i— 1 Cycle 

where fi and pi are the position and momentum of particle i respectively. The only parameter of this Hamiltonian is 
the energy density of a fundamental flux tube. Using the dimensionless variables introduced in Sec. IVIl Eq. (123)) can 
be recast in the form 

4 

r H = ^2^Jq^+ ^2 \xi-x i+ i\, (26) 

i—1 Cycle 



FIG. 3: Schematic representation of the four-gluon system considered. Four fundamental flux tubes (gray lines) link the 
constituent gluons (circles), illustrating Hamiltonian (|26|) . The relative coordinates defined by relations (|28[1 are also represented 
(black arrows), while the center of mass, denoted as R, is given by Eq. (|27[1 . The (1,2) and (3,4) gluon clusters have been 
surrounded. 



in which no free parameter is present. 

Hamiltonian (J2SJ) seems too simple at first sight. Indeed, pair interactions coming from one-gluon-exchange processes 
have not been included. However, it has be shown in Ref. [37| that the relativistic corrections to bring to the confining 
potential are proportional to two color factors. One of them is a projector on a color octet state and vanishes in our case 
because a gluon pair is always in a color (anti)decuplet. The other color factor also vanishes because it is proportional 
to C gg — 6, which vanishes in the color (anti)decuplet channel. Remarkably, the few physical ingredients that have 
been put in Hamiltonian (|26|) appear to be enough to accurately describe the H glueball thanks to the vanishing 
of the relativistic corrections in this case. It is possible that other nonperturbative effects like instanton-induced 
interactions could contribute in a non trivial way, but we neglect them in this simple first approach. 

For future computations, it will be very convenient to have an appropriate set of relative coordinates at our disposal. 
We first define R the center of mass of the system as 



R 



A Z\2 Xi > 



and the three relative coordinates {X, Y, Z} as 



i=l 



(27) 



X = Xi~X2, Y = X3~Xi, 



z 



Xi + x 2 




2?3 + X4 


2 




2 



(28a) 
(28b) 



The meaning of X and Y is obvious, while Z is the separation between the center of mass of the cluster made by 
gluons 1 and 2 and the corresponding one for gluons 3 and 4. A graphical representation of the four-gluon system 
and the relative coordinates is given in Fig. [3] 

We have shown in Sec. lIV Cl that the internal structure of the H can be elucidated thanks to the helicity formalism 
provided that the particular permutation symmetry of this state is taken into account. A crucial point has now to 
be outlined: By definition of the relative variables X and Y , the orbital angular momenta corresponding to clusters 
(1,2) and (3,4) are nothing else than Lx and Ly respectively. Consequently, we are led to the following results 



L 



x 



(29) 



These highly non standard values for a square orbital angular momentum are a direct consequence of the helicity 
formalism and could never have been deduced from an usual, spin-based, approach. A look at definition (|28b|) then 
shows that Lz can be identified with the relative orbital angular momentum of the two clusters, so we obtain 



f 2 
^ z 



0. 



(30) 



B. Numerical results 



We have now all the necessary ingredients to compute at least a crude approximation of the mass of the lightest 
H state within our model. The mass that we will obtain can be compared to the lattice QCD data 

r M +- = 11.54 ±0.24 0, r M +- = 11.66 ± 0.19 0. (31) 

Both results are actually compatible, the second one coming from a more recent computation. 

The glueball mass we look for should a priori be numerically computed from Hamiltonian (|26p . A full four-body 
calculation with such a semirelativistic Hamiltonian is a very difficult task in itself, but here the situation is even worse 
because the helicity formalism imposes non standard values for the matrix elements of the orbital angular momentum. 
Since the four-body computations are widely out of the current numerical methods at our disposal, we rather choose 
a "digluon-digluon" approach: We will first compute the mass of the (1,2) and (3,4) clusters, and then calculate the 
mass of the bound state made of these two clusters. This scheme is very useful because only two-body calculations 
are needed, but also because we have investigated the internal structure of the +_ state in this particular scheme 
throughout this paper. Both gluon clusters should actually have the same mass ro M c , given by the lowest-lying state 
of Hamiltonian 

r H c = 2J$%+X, (32) 



where qx and X are conjugate variables and where X = \X\, with \J->x) = 4 f rom Eq. (|29[) . Such a Hamiltonian can 

be readily understood from Fig. [31 Equation ([3"2"|l actually refers to the (1,2) cluster, but Hamiltonian 2-y/^T + Y, 
valid for the (3,4) cluster, would lead to the same mass by symmetry. The eigenvalues of Hamiltonian (132]) can be 
accurately computed thanks to the Lagrange mesh method [39(; we find ro M c = 4.72. 

As the masses of the clusters are known, the total mass ro M of the system can be computed. It should be given 
by the ground state of Hamiltonian 



r H g = 2y/% + rlMl + 2Z. (33) 

One can indeed conclude from Fig. [3] that the two "digluons" should feel a confining potential such as if they were 
bound by two fundamental flux tubes. Applying again the Lagrange mesh method for {^z^ = [ see Eq. ([3H]) ]. one 
has 

r M = 11.61, (34) 

in nice agreement with the lattice QCD data (|31| . 



VIII. CONCLUSIONS 



In the present work, we have given arguments showing that the glueballs and gluelumps that are currently observed 
in lattice QCD can be understood in terms of bound states of a few massless constituent gluons with helicity- 1. In 
this scheme, the lowest-lying C = + glueballs can be identified with two-gluon states (or at least with hadrons in 
which a two-gluon component widely dominates) , while the lightest C — — glueballs are compatible with three-gluon 
states and the heavy H glueball should be seen as a light four-gluon state. For what concerns gluelumps, we have 
shown that the lightest ones, which are C = — states, are compatible with one-gluon states. This is allowed because 
of the static color octet source. To confirm this classification, we have shown that a linear behavior of the masses with 
respect to the expected number of constituent gluons is observed, suggesting that each constituent gluon brings a 
given energy to the total mass. The slope of the mass formula can be reproduced within a simple flux tube model with 
linear confinement and one-gluon-exchange potential provided that the gluons helicity is correctly taken into account. 
A useful concept has been introduced in order to classify the various gluonic states, that is the natural parity, defined 
by (— l) Ns+J . It appears that a glueball with natural parity will always be lighter than the corresponding glueball 
with the same J and C but the opposite parity. 

It is worth recalling that, if constituent gluons were regarded as spin-1 particles, some states would appear that are 
not observed in lattice QCD. The most striking ones would be light 1 ±+ glueballs in the two-gluon sector, but also a 
gluclump and a light three-gluon glueball. The absence of these states clearly favors helicity-1 constituent 
gluons. 



In order to convince ourselves that the 0^ glueball is a relevant four-gluon candidate, we have computed its mass 
starting from a four-body flux tube Hamiltonian. Due to the unusual color wave function and permutation symmetry 
of this glueball, only the confining potential is present and other relativistic corrections vanish. Other mechanism as 
instanton-induced interactions have been neglected. Moreover, the assumption that the constituent gluons are helicity- 
1 particles leads to the knowledge of the internal structure of the considered glueball. Once all these ingredients are 
put together, the mass we obtain is in agreement with lattice QCD without fitting any parameter. This brings a good 
support to the constituent-gluon picture we developed in this work. 

The above results are actually a strong motivation to reconsider the various potential models of hadrons with a 
gluonic content under the assumption that constituent gluons are massless particles with a helicity-1. If two-gluon 
glueballs have already been studied in Ref. [ill ], the mass spectra of gluelumps, hybrid mesons, and three-gluon 
glueballs should now be explicitly computed within the full helicity formalism. Such a work is in progress. 
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APPENDIX A: HELICITY STATES FOR GLUELUMPS 

Using the formalism developed in Ref. [HI], it is easy to build the helicity states for a gluelump, that is a bound 
state with a gluon and a static ++ color octet source, with J > 1, 

|1,0;J P ,M,1) = |3j_i j) + ^_jL_ |3 J+lj) with p =(_)', (Al) 

|1,0;J P ,M,-1) = -| 3 Jj) with P = (-) J + 1 . (A2) 
With these decompositions, a simple computation gives (L 2 ) = J (J + 1). 



APPENDIX B: C-PARITY OF FEW-GLUON GLUEBALLS 

The gluonic field is defined as A p = A p T a , where T a = A a /2 are the generators of su(3) and where A Q are the 
Gell-Mann matrices. If C is the C-parity operator, the charge conjugate of A^, that we write A M , is such that [Hj] 

A ii = CA il C- x = -A T ll . (Bl) 

The index T denotes a transposition. Let us begin with the simple case of two-gluon states. The only way to make 
a color singlet is to write such a state as 5 a bA a ^A b u , S a b being the Kronecker symbol. As Tr(A a A/,) = 25 a b [Hj]> we 
have 5 a t,A®A°, oc Tr(A M j4„). It is then readily checked that Tr(A fl A v ) — ^(A^A^), implying that two-gluon states 
always have a positive C-parity. Next, there are two ways of making a color singlet with three gluons: fabcA^A^Ap cx 
Tr([A M , A v \ Ap) and d a bcA a pA b v A c p cx Tt({A^,A v }A p ), with f a b c the totally antisymmetric structure constant of su(3) 

and d a bc the totally symmetric ones. These are actually the [[8, 8] 8a , 8] la and [[8, 8] 8s , 8] ls configurations respectively. 
It can then be computed that 

Tr([i M , A„]A P ) = Tr([j4 M , A v ] A p ) => C = +, (B2) 

Tr({i M , A V }A P ) = -Tr({^, A v }A p ) =► C = -. (B3) 

These configurations have the permutation symmetries [i^vp] and {[ivp) respectively. 

The case of four gluons is a little bit more involved. Let us first focus on the particular color wave func- 
tions in which the gluon pairs are in a color octet. The following combinations lead to a color singlet: 
fabcf c de A-A b v A d p A% cx TiHAp, A u ][A p ,A a }), d abc d c de A-AlA d p A% cx |lV({{A M , A v }, A p }A a ) - iTr^A^Tr^AO, 



and f a } lc d c de A a ^A b 1J A d p A e a oc TrQA^, A„]{j4 p , A a }). They correspond respectively to the color wave functions 
[[8,8] 8 ", [S,^ 8 *] 1 , [[8,8] 8 -, [8,8] 8s ]\ and [[8, 8] 8s , [8, 8] 83 ] 1 . As previously, the C-parity can be computed. First, 

TrQA^AuliA^Aa]) = Tt([A p , A U ][A P , A a }) => C = +, (B4) 

12 21 

with a permutation symmetry that can be of the form Mp^pa only, the total antisymmetry [pvpcr] being excluded by 
virtue of the Jacobi identity. Second, we have 

^Yi{{{Ap,A u },A p }A a ) ~ ^{A^AjTriApAv) = ^({{^, A v ], A p }A„) - ^Tv(A^)Tv(A p A a ) 

=> C = +, (B5) 

12 21 

with a permutation symmetry that can be of the form (pi/ pa) or Mpv pa ■ The only configuration that can lead to a 
negative C-parity actually is 

TT([A p ,A u ]{A p ,A a }) = -Tr([A^ A„]{A P , A a }) => C = -. (B6) 
A priori, both the permutation symmetries Mpfu pa and M^Jpa are allowed, but the equality f a ^ b c d de j c = forbids 

13 ll 

the Mpvpa symmetry. The last equality just reflects the obvious fact that, for any three given matrices A , B and C , 
the following identity holds [A, {B, C}] + [B, {A, C}] + [C, {A, B}} = 0. 

The two remaining four-gluon color wave functions are now [[8, 8] 10 ° , [8, 8] 10 "] 1 and [[8, 8] 27s , [8, 8] 273 ] 1 . The study 
of these configurations can be achieved by building the various representations we need from bilinear quantities in the 
fundamental irreducible representation of su(3) and its conjugated, 3 and 3 , using 3 ® 3 = 8 © 1 [43(. Let us take a 
complex 3-vector Vj and its complex conjugated V k transforming in the 3 and 3 representations respectively, where 
j, k = 1,2,3, and impose the vanishing of the trace V l Vi = . Then, T\ = i V^Vk defines an anti-hermitian traceless 
3x3 matrix, therefore transforming in the adjoint irrep of su(3), (T^)* = -i {V°Vk)* = -i {V k V.j) = -T k . This 
relation is actually related to Eq. (IB1|) . Note that the dual irreps of su(N) are equivalent to the complex conjugated 
ones, implying for su(3) that 

U? A , m ~ Ai B s e^ k (B7) 

belongs to the 3 representation, the symbol ~ meaning 'isomorphic to'. The tensor product 8 ® 8 can then be 
performed explicitly, giving the well-known result 

8<g>8 = 27 s ©T0 a ®10 a ®8 a ®8 s ©l s , (B8) 

where the permutation symmetry, i.e. symmetric (s) or antisymmetric (a), of all the representations have explicitly 
been written. Starting from (A l Aj)(B k Bi), the left-hand-side of the above relation, we find the correct decomposition 
with 

l s ~ Z(A, B) ee (A'AgX&Bp) , 8 S ~ Y)(A, B) = (A i A p )(B p B j ) + {&B v ){APAj) - | 5) Z(A, B) , 
8 a ~ (A i A k )(B k Bj) - (B^B^A,) , 10 a ~ A«B k U$ B] , 10 o ~ U^AjB^ , 
27 s ~ X$(A, B) ee X$(A, B) - 1 *grj> (A, B) - 1 5$6*j Z(A, B) , 

with Xjf(A, B) — A^B k ^AfjB^ and B) = 0. The action of the C-parity on all these irreducible representations 

is particularly simple: It lowers (raises) the upper (lower) indices and removes (adds) bars accordingly. For example, 
Aj is changed into AK It is readily checked that the l s has a positive C-parity, and that all the previous results 
concerning four gluons are recovered. In the case of a four-gluon system, the decomposition (|B8|) tells us that two 
singlets can be obtained for gluon pairs in 10 a and 10 Q . These are given by 

1 ± * U$ B] A jBl) u{™ D] &D» ± U^m^C jDl) . (B9) 

This last line can be used to check that the 1+ (1_) singlet has a positive (negative) C-parity. Finally, their 
permutation symmetry is of M' 2 ' 2 ' type. The configurations involving the 27 s representation generate repulsive 
forces between the gluons. They will necessarily lead to heavy glueballs and will consequently not be studied here, 
were attention is drawn on low-lying states. 



A few comments can be done concerning gluelumps. The color octet source can actually be seen as a color scalar 
particle (f> — (f> a T a such that (j> = (fi T . Then, a one-gluon gluelump corresponds to Tr^A^), and it is easily checked 
that the C-parity of such a state is always negative. For two-gluon gluelumps, one can further compute that C = + 
(— ) for a totally symmetrical (antisymmetrical) color configuration. 
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